in 
o 
o 

(N 



43 

a: 



- 1—1 

X 



Complete CKM quark mixing via Dimensional Deconstruction 



P. Q. Hun£ 

Institute for Nuclear and Particle Physics, 
and Department of Physics, 
University of Virginia, 
Charlottesville, Virginia 22904-4714, USA 

A. SoddJ] 

Department of Physics, 
National Taiwan University, 
CS| ! Taipei, Taiwan 106, R.O.C. 

m : n 

! Ngoc-Khanh Trar|i| 

ON 



Department of Physics, University of Virginia, 
Charlottesville, Virginia 22904-4714, USA 



q i (Dated: February 2, 2008) 



Abstract 



^h. It is shown that the deconstruction of [SU (2) x U (V)] N into [SU(2) x U(l)] is capable of providing 



all necessary ingredients to completely implement the complex CKM mixing of quark flavors. The 
hierarchical structure of quark masses originates from the difference in the deconstructed chiral 



zero-mode distributions in theory space, while the CP-violating phase comes from the genuinely 
complex vacuum expectation value of link fields. The mixing is constructed in a specific model to 
satisfy experimental bounds on quarks' masses and CP violation. 
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I. INTRODUCTION 

Dimensional deconstruction 

US 

is a very interesting approach to dynamically gener- 
ate the effects of extra dimensions departing from the four-dimensional (4D) renormalizable 
physics at ultraviolet scale. That is, apart from having the viability in the sense of renor- 
malizability, whatever amusing mechanisms being dynamically raised by the virtue of extra 
dimensions (ED) now can also be easily arranged to rise dynamically in a pure 4D framework. 
In this paper we look specifically into two such important mechanisms of extra dimension 
theories, namely the localization of matter fields in the bulk U, 0] and the dynamical 
breaking of CP symmetry by ED Wilson line □ □ □ . Ultimately, the hybrid of these 
two mechanisms is just the well-known complex mixing of fermion flavors. And it is concep- 
tually interesting to note that dimensional deconstruction (DD) nicely encompasses both of 
these issues. In other words, complete Cabibbo-Kobayashi-Maskawa (CKM) mixing can be 
generated naturally via dimensional deconstruction. 

With the presence of extra dimensions, one has a new room to localize the matter fields 
differently along the transverse directions as in the so-called split fermion scenario. Various 
overlaps of fermions of different flavors then induce various fermion masses observed in 
nature (see e.g. Amazingly, the deconstruction interaction is also able to 

produce similar localization effects Indeed, after the spontaneous breaking of link fields, 
fermions get an extra contribution to their masses via the Higgs mechanism. Fermions 
then reorganize themselves into mass sequences and the lightest mass eigenstate of these 
towers exposes some interesting "localization" pattern in the theory space (also referred to 
as deconstruction group index space). We will first work out the analytical expressions and 
confirm the localization of these zero modes in a rather generic deconstruction set-up. The 
next question to raise is how to make these light modes chiral. Imposing some kind of chiral 
boundary conditions [2| is the answer again coming from the ED lessons. There is however 
one more subtle point to be mentioned here. If one truly wishes to relate the ED scenario 
to the dimensional deconstruction, one needs to latticize the extra dimensions to host the 
deconstruction group. There comes the lattice theory's issue of fermion doubling, and its 
standard remedy, such as adding to the Lagrangian a Wilson term 1J| would remove half of 
original chiral degrees of freedom. This is the reason why most of previous works addressing 
the fermionic mixing in deconstructed picture (e.g. [3,13,13) usually start out with only 
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Weyl spinors. In the current work, we adopt a different and somewhat more general 4D 
deconstruction approach where no extra dimension is actually invoked. As a result the 
fermions to begin with keep a standard 4-component Dirac spinor representation. 

In any deconstruction set-up, the link fields transform non-trivially under at least two 
different gauge groups. This implies a complex vacuum expectation value (VEV) for these 
fields, whose phase would not be rotated away in general. After the deconstruction pro- 
cess, this phase is carried over into the complex value of wave functions and wave function 
overlaps of fermions. In turn, the induced complex-valued mass matrices can render a re- 
quired CP- violating phase in the well-known KM mechanism. In contrast, we note that 
the generation of complex mass matrices within the split fermion scenario is a non-trivial 
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problem and requires rather sophisticated techniques to solve 
above CP violation induction via deconstruction can also be visualized in extra dimensional 
view point. Indeed, because of having the same symmetry transformation property, DD link 
field can be identified with the Wilson line pointing along a latticized transverse direction 
(Appendix B), and the latter then can naturally acquire a complex VEV in the generalized 
Hosotani's mechanism (f| 3, S of dynamical symmetry breaking. Apparently, the source 
of CP violation in this approach comes from the complex effective Yukawa couplings so it 
can be classified as hard CP violation. Nevertheless, those couplings acquire complex values 
after the spontaneous breaking of the DD link fields. In that sense this CP violation pattern 
could also be considered soft and dynamical. 

This paper is presented in the following order. In Sec. IIA we give the zero mass eigenstate 
of fermions obtained in the deconstructed picture, in Sec. IIB the resulting expression of 
mass matrix elements, and in Sec. IIC the symmetry breaking of [577(2) x U(1)] N into 
[SU(2) x Z7(l)]. In Sec. Ill we present the numerical fit for quark mass spectrum and CKM 
matrix in a model where each "standard model" Higgs field is chosen to transform under 
only a single deconstruction subgroup. The conclusion and comments on numerical results 
is given in Sec. V. Appendix A provides a detailed derivation of zero mode wave functions in 
4D deconstruction using combinatoric techniques. Appendix B outlines intuitive arguments 
on the complexity of link field inspired by lattice models. Appendix C presents analytical 
expressions for wave function overlaps used in the determination of mass matrix elements. 
Finally, Appendix D gives referencing values of key physical quantities that have been used 
in the search algorithm (Table I), and numerical solution of our models' parameters (Table 
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II). 



II. DECONSTRUCTION AND QUARK MASS MATRIX 

In this section we describe how the mixing of quark flavors arises in the DD picture. But 
we first briefly recall the basic idea of the dimensional deconstruction applied to just a single 
quark generation. The family replication will be restored in the later sections. 



A. Zero-mode fermion 



We begin with N copies of gauge group [SU(2) x U(l)] n where n = 1,..,N. To each 
group [577(2) x U(l)] n we associate a 5?7(2) n -doublet Q n , and two 5?7(2) n -singlets U n , D n . 
These fields transform non-trivially only under their corresponding group [577(2) x U(l)] n 
as (2, qo), (1, qu), (1, qn) respectively, with g's denoting ?7(l)-charges. Finally, we use 3(7V — 
1) scalars 4>n-i,w 4>n-i,n, 0n-i,n transforming respectively as (2, q Q \2, -q Q ), (1, qu\l, -qu), 
(l,q D \l,-q D ) under [SU(2) x C/(l)] n _i x [SU(2) x U(l)] n to "link" fermions of the same type. 
Because of this, scalars 0's are also referred to as link fields hereafter. For the simplicity of 
the model, we assume a symmetry for the Lagrangian under the permutation of group index 
n. 

The n^Li [577(2) x U(l)\ n gauge-invariant Lagrangian of the fermionic sector is 

/ N N-l N \ 

C = {J2 QniTPnQn + £ Qn<f>in + lQn+l ~M Q J2 QnQn + {Q U) + (Q D) , (1) 
\n=l n=l n=l / 

where n denotes the covariant derivative associated with gauge group [577(2) x C7(l)] n , 
and Mq, Mjj, are the bare masses of fermions. Ultimately, we are interested in achieving 
chiral fermions of standard model (SM) at low energy scale. To this aim we impose the 
following chiral boundary conditions (CBC) on fermion fields Q]. 

Qir — Qnr = , 4>n-i,nQn,l = VqQn-i,l , 

U\l = Unl = , (/>n-i,nUn,r = VuUn-i,r, (2) 

D\l = Dnl = , n-Dn,r = VdDn-i,r ■ 



We note that those conditions are in agreement with the gauge transformation property of 
fields, e.g. (fi%_ ltN QN,L and Qn-i,l transform identically under the underlying gauge groups. 



Essentially, these boundary conditions render one more left-handed degree of freedom over 
the right-handed for Q field, and the contrary holds for U and D fields. The actual cal- 
culation will show that the zero-mode of Q field indeed is left-handed while for U, D it is 
right-handed. When the link fields <fiQ> u ' D assume VEV proportional to Vq ud, above CBC 
become the very reminiscence of Neumann and Dirichlet boundary conditions. 

In the deconstruction scenario, after the spontaneous symmetry breaking (SSB) the link 
fields acquire an uniform VEV Vqjj,d respectively, independent of site index n (in accordance 
with the assumed permutation symmetry), and the fermions obtain new mass structure. 
Using the CBC (J2J), the fermion mass term can be written in the chiral basis as 



I QlL ^ 



(Q 



2R, 



Q 



N-l,R, 



+ (Qil,..,Qn-i,l)[M q ]* 



\ Qn-x,l j 
(Qr,l <-> U L>R ) + {Qr, l 



Q2R 
\ Qn-i,r J 



D 



L,R) j 



(3) 



where the matrix [Mq] of dimension (N — 2) x (N — 1) is 

f -V% M Q -V q 



\(N-2)x(N-l) 



-Vq Mq -Vq 
-V% Mq 



\ 



Mq - 
-Vq* M q 



■Vo 



V, 



Q j 



(4) 



By interchanging Qr^l Ul,r, Qr,l Dl,r, the matrices [Mu], [Md] of dimension (N — 
1) x (N — 2) can be analogously found. 

By coupling the following Dirac equations for chiral fermion sets {Qr} = 
(Q2R,--,Qn~i,r) T and {Q L } = (Q 1L , .., Q N - 1:L ) T 



^{Qr}-[M q ]{Q l } = 0. 



i0{Q l }-[MqV{Q r } = O. 



(5) 



we see that [MqMq] is the squared-mass matrix for the left-handed components Ql and 
[MqMq] for the right-handed Q R . Since at low energy, we are interested only in the chiral 
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zero modes of fermions, we will work only with [MqMq], [MjjM^], [MdMjj] in what follows. 



( Wq? 
-m q v* 





[M+Mq] 



-MqVq 
M* Q + 

W Q \ 2 



-2M a V, 



Q 



(VA) 2 



V 2 Q 
-2M Q V Q 

2|Vq| 2 
-2M n VA 



\ 



-2M Q V Q 

M Q + 
2\Vq\ 2 



-2M Q V* 
(V*) 2 



M Q + 
2|Vq| 2 



-2MqVq 



-2M Q V Q 

M Q + 
2|Vq| 2 



(VA) 2 



V2 Q 

-2M q Vq- 
Wq\ 2 

-2M Q V*+ M*+2\Vq\ 



(Vq + V*)M q J 



(6) 

and similar expressions hold for [MjjM}j] , [MdM^] . The quantitative derivation of the zero- 
eigenstates, which are identified with the SM chiral fermion, is presented in appendix A. In 
this section we just concentrate on some qualitative discussion. In general the diagonaliza- 
tion of matrices (JBJ) leads to the transformation between gauge eigenstates Q n i and mass 
eigenstates Q m L 

QnL = %^Q\nmQmL i QnL = [^olmnQmL , (7) 

where the matrix \Uq] diagonalizes [MqMq] 

W ] Q M Q ] diag = \Uq\ * [MqMq] [Uq] . (8) 

The key observation, which will be analyzed in more details in appendix B, is that VEV 
Vq,u,d are generically complex and \Uq^u,v] are truly unitary (i.e. not just orthogonal). This 
in turn gives non-trivial phases to zero-mode fermions Qol,Uql, Dql in Eq. ((7j) and after 
the SM spontaneous symmetry breaking the obtained mass matrices are complex. Further, 
the explicit solution of zero mode Qql (and Uor, Dor) in the mass eigenbasis exhibits a very 
interesting "localization" pattern in the group index space n (see Appendix A). This in turn 
can serve to generate the mass hierarchy among fermion families in a manner similar to that 
of ED s pUt fami0n — (see ,, QQ). T tas we S ee that di—a, de—tion 
indeed provides all necessary ingredients to construct a complete (complex) CKM structure 
of fermion family mixing. 



B. Complex Mass Matrix 



In order to give mass to the above chiral zero-mode of fermions, we introduce Higgs 
doublet fields just as in the SM. In the simplest and most evident scenario (see ^3), there 
is one doublet Higgs H n transforming as (2, qg — q D = q v — qg) under each \SU{2) x 
U(l)] n group. We also implement the replication of families by incorporating family indices 
i,j = 1,..,3. Another scenario to generate the (vector-like) fermion mass hierarchy by 
assuming various link fields to connect arbitrary sites of the latticized fifth dimension has 
been proposed in l^ . 

The gauge-invariant Yukawa terms read 



N 



N 



4 £ Qfi^KU^ + 4 £ Q®H n Djp + H.c. . 



(9) 



n=l 



71=1 



In order to extract the terms involving zero modes, which are the only terms relevant at 
low energy limit, we rewrite Q in the mass eigenbasis. However, this procedure depends 
explicitly on the specific CBCs being imposed on each of the fields Q, U, D. To be generic, 
let us consider the following configuration. We assume the "localization" of zero modes Qol, 
Uqr and Dqr to be at n — 1, n = 1 and n = N respectively. To achieve this localization 
pattern, we impose the following CBCs on these fields (see Eq. (J2J) and Appendix A, Eq. 

21) 



(10) 



Because of these boundary conditions, zero modes Qol, U or and D 0R would be localized at 
n — 1, n — 1 and n = N respectively, this also means that the first term of Eq. © would 
represent the overlap between 2 wave function localized at the same site n = 1, while the 
second term represents the overlap between wave functions localized at n — 1 and n = N . 
Using (|10|) to eliminate the dependent components and after the SM spontaneous symmetry 
breaking (H n ) = (0,v/^/2) T uniformly for all n's, we can rewrite the Yukawa term (jOJ) as 



QlR. 


- o (i) 


= o, 


tN—1,N^CNL 


— V Q WN-1L 


rrU) 
U 1L 


- u (j) 


= 0, 


( PN-1,N u NR 


— v v u N _ 1R , 


D (k) 


- D {k) 

— U NL 


= 0, 


,(fc)Dt (fc) _ 
( Pl,2 U \R ~ 


- i/( fc )* n( fc ) 

" V D LJ 2R • 



m=i 

3 

V 2. i k=l 



.D 



V 



r(k)* 



U) 



^CN-1L U N-IR + WnL u nR, + { ^nR u nL 



r(7) 



Af-1 



di)Q Aj)U 
'N-l,N tN—1,N 



W ttU) 



D ' ^(i) n (fc) , 
' K 4\L U 2R ' 



Vr 



n=2 
Af-1 



+ 



h,2 



di)Q 
J N~1,N 



^N-lL^NR + WnL^nR ' ^nR^ 



nL) 



n=2 



After going to the mass eigenbasis by the virtue of transformation of the type (jjj), keeping 
only zero-mode terms and together with the assumption of universality for the Yukawa 
couplings in the up and down sectors, we obtain the following effective mass terms 



E + e of L M? k m , 

i,j=l i,k=l 



(12) 



with 



V2 



fN-2 



nS&U \n,0 + 



V, 



V, 



u 



vn=l 



iMQ AW 

J N-1,N YN-l,N, 



(i) 

Q IN 



-i,o[^r/ W-i,o 



(13) 



(&)* 



>1,2 



Mi 



+ ( E [^ ) ]; [^]n,0 ) + Wl^o I + "^g— P< 



/(*) 



(0i* 



r( k h 



V, 



(0 



i n=3 



)Q 

N—1,N 



(14) 



Because all [ZYq], [ZYy], [Wd] are unitary the mass matrices M", M d are generally complex. 
Thus in this simplest deconstruction approach, we might better understand the dynamical 
origin of CP- violation phase in the SM mass matrices. We also note that (|13p. (|14|) represent 
the specific case where Qol, Uqr and Dqr are localized at n = 1, n = 1 and n = N 
respectively. All other localization configurations can be similarly found. Further, when we 
replace link fields 0's in ()13j) . (|14jl by their VEVs following the deconstruction, these mass 
matrix elements will look much simpler (see ([27)1 . (|28p). 

Before moving on to give explicit expressions of these complex-valued mass matrices in 
term of zero mode wave functions (appendix A) and perform the numerical fit, let us briefly 
turn to the breaking pattern of product group Yl^ =1 [SU(2) x U(l)] n . 



C. Deconstructing [SU{2) x U(1)] N 

For the sake of completeness, in this section we will describe the breaking of [577(2) x 
U(1)] N into the SM [SU(2) x U(l)} gauge group by giving uniform VEVs to link fields. 
The transformation and charge structure of fermions and scalar link fields have been defined 
in the beginning of previous section. To identify the unbroken symmetries following the 
deconstruction, we look at the covariant derivative and kinetic terms of scalars. 

D »<t>n,n+1 = ^n,n+l ~ i( lU y ^n/^n.n+l + ^l/y B n+lfi(p n ,n+l 7 ( 15 ) 
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D, 



kD 



g'o 



9o 



(16) 



Mr n,n+l "^rn^ta - 2 ""f* T ti,t\-\-x ' - ^ 

where B n is the gauge boson associated with U{l) n) while g' is the common gauge coupling 
for all [/(l)'s. For abelian groups, the opposite signs of the last two terms in (fi3j) (and also 
in (JIHl)) originate from the opposite charges of </>^ n+1 ( an d 4>nn+i) under U(l) n and U(l) n+ i 
(so that terms like U n (p^ n+1 U n are gauge-invariant). 

For non-abelian groups, the similar sign reversing will hold for terms in the expression 
of covariant derivatives (see Eq. (}2T]) below), the nature of which also has its root in the 
gauge invariance of the theory. Indeed, under the Yang-Mills SU{2) n x SU {2) n+ i gauge 
tranfromation (note that (j)n, n +i is a 2 x 2-matrix) 



T 

- 1 n. 



n,ri+l ' - L nH J n,n+l^ n +l i 

Qn * T n Q n , Q n +\ 
T ... r T 



.4 



rt/i , 



T 



A 



nfi 



A 



n+l/i ~ 



T 



n+1 



r 



n+i/i ; 



2"n+l ( , 9/i^n+l)7 1 7i+l 
#0 



(17) 
(18) 
(19) 

(20) 



the covariant derivative of 4>nn+i m ust be formulated as follows (so that it transforms exactly 



like 0^ n+1 in (HID) 



Dfi4>n,n+1 — <9/J07i,n+l (^Q ^^7i/j0n,n+l + 'SoAi/j^Rjn+b 



So 



r , 



+(^Qy 5 n+i,,0„,77+i + iflb^ft+iAvj 



(21) 



where A n and T n are respectively the gauge bosons and some 2 x 2-special unitary matrix 
characterizing the SU(2) n transformation, while go is the common gauge coupling for all 
S77(2)'s. 



After the deconstruction 



lU,D 
'n,n+l 



U,D, V n ,n+1 



V, 



Q 



L 2X 2, the mass terms for gauge 
bosons are generated. Specifically, we obtain as parts of kinetic terms (D /J 0^ n+1 )^(D /i 0^ n+1 ), 
(D^Zn+iVi^^n+i)^ ^A( D ^4>n,n+i)K D ^4>n,n+i)} the following gauge bosons squared mass 
matrices 



/' 1 -1 
-1 2 



[M%] = A 



B 



\ 



2 -1 
-1 1 J 



( 1 -1 
-1 2 



\ 



2 -1 
-1 1 J 



(22) 



where, after restoring the family replication index (i = 1, 2, 3), 

Ab = E^I^ ) I 2 + ^I^ ) I 2 + ^I^ ) I 2 ), A^Esol^fl 2 - (23) 

1 1 

Both matrices in (}2*2*j) have a "flat" zero eigenstate. This indeed indicates the uniform 

breaking of [SU(2) x U(1)] N into the diagonal (SM) group [SU(2) x U(l)], whose gauge 
bosons are massless and given by 

1 N - i N 

5 M = -7jtf£A**> 4 = ^E^- (24) 

v A n=1 VJ\ n=l 

In Eq. ()23|) it is also shown that the pattern of symmetry breaking is not spoiled by family 
replication as long as charges qy (and qo, qQ ) are independent of the site index n under a 
presumed permutation symmetry (just like V^%q). Finally, by extracting the interaction 
between fermions and massless gauge bosons from fermion kinetic terms in JI} one can see 
that the couplings of the unbroken group scale as g' = g' /VN and g = go/\H\f, while the 
charge structure (of fermions in mass eigenbasis) under this diagonal group remains intact. 



III. FITTING THE MODEL'S PARAMETERS 
A. Model, Parameters and Numerical Method 

In the previous section we have outlined the process diagonalizing the squared-mass 
matrix (JHJ). The complete diagonalization process is complicated, but as we are concerned 
only with the zero eigenvalue problem, the computation can be done analytically in the 
general term (see Appendix A). Since [Uq] diagonalizes [MqMq] ©, the zero eigenstate of 
[MqMq] is just the first column of \Uq], i.e. in the notation of appendix A 

[Wq ] n ,0 = Xq u , (25) 

and similarly 

[Z^ko = xf n , [U^U = y { nl , (26) 

where x n 's are given in (|A20|) (corresponding to a zero mode localized at the end point 
n — 1) and y n 's in (|A24J) (corresponding to a zero mode localized at the end point n = N). 

After the spontaneous symmetry breaking, the link fields acquire an uniform VEV Vq^d 
respectively (independent of site index n). In term of x n 's and y n 's, the SM mass matrices 
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(JIIIJ), (HU) for up and down quark sectors become 




v lYv^ 1 W* C?A i (*)* 



(3) 



(27) 





(28) 



where sc n 's, y^s are given in (|A20|) . (|A24|) respectively. The analytical forms of Ij27j) . (j2Bj) in 
term of model's parameters are worked out in Appendix C, Eqs. ()C1|) . (|C2|) . 

Again, let us remind ourselves that (J2*7)l represents the overlap between two wave functions 
localized at the same site n = 1 while (|28|) represents the overlap between one wave function 
localized at n = 1 and the other at n = N. The model under consideration consists of 20 
real parameters (see Table II and appendix B): 3 complex VEV V's for each complete quark 
generation (Q,U,D)i (i=l,2 or 3), and 2 real "dimensionful" Yukawa couplings Kuv/y/2, 
We choose to fix = 10 throughout. 

The numerical approach to fit the parameters consists in minimizing a positive function 
which gets a zero value when all the predicted quantities are in the corresponding experimen- 
tal ranges Q|. The minimization procedure is based on the simulated annealing method, 
which seems workin g be tter than other minimization approaches when the parameter space 
becomes larger l^.^ll]. The input referencing physical quantities are given in Table I of 
Appendix D. 

We consider eight different cases, which correspond to all the eight possible ways of 
localizing the left and right components. The eight different cases are the following: 

1. Q, U and D localized in n = 1 denoted as (QUD1) 

2. Q and U localized in n = 1, D localized in n = N denoted as (QU1DN) 

3. Q and D localized in n = 1, U localized in n = N denoted as (QD1UN) 

4. Q localized in n — 1, U and D localized in n = N denoted as (Q1UDN) 

5. Q, U and D localized in n = N denoted as (QUDN) 

6. D localized in n = 1, Q and U localized in n = N denoted as (D1QUN) 

7. U localized in n — 1, Q and D localized in n = N denoted as (U1QDN) 
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8. U and D localized in n = 1, Q localized in n = N denoted as (UD1QN) 

We specially note that, due to the mirror complexity between CBCs (0) and (jA22|) . the 
mass matrices obtained in the cases 1 and 5, cases 2 and 6, cases 3 and 7, cases 4 and 8, are 
complex conjugate pairwise. In the result, all eight cases are inequivalent. 



B. Numerical results 

In the following we present the characteristically important numerical results for the four 
cases out the eight mentioned above, for which we were able to find solutions. The cases are 
referred to in the above order. For each case we give one particular, but typical, numerical 
complete set of the 20 defining parameters (Table II), the quark mass matrices and quark 
mass spectra, the CKM matrix and the CP parameters. Complex phases are measured in 
radiant, and N = 10 for all cases. The masses are given in GeV and are evaluated at the Mz 
scale. For the sake of visualization, we also present graphically the comprehensive solutions 
of the quark wave function profiles in the theory space (Fig. the mass spectrum (Fig. 
12} , the CKM matrix (Fig. EJ) and the p-f] CP parameters (Fig. for the case of all fields 
Q, U and D localized at the same site n = 1. 

• Case 1: Q, U and D localized in n = 1. 



M (Qum) = 78AGe y 



m (QUDi) = _ 0021 GeV _ 



0.925 e 



-0.558? 



0.027 e 2m9i 

,0.306? 



\ 0.948 & 



0.923 e~°- 501i 
0.027 e 2Mei 
0.942 e 0367i 



0.951 e 



-0.570i \ 



0.029 e 2 006i 
0.973 e a280i 



m 



(qudi) 



0.702 GeV, m\ 



(QUDI) 



181.1 GeF. 



(29) 



(30) 



M, 



(QUDI) 



,1 



0.909 e 



1.490?' 



0.782 e 



-2.072? 



0.960 e 



0.992? \ 



V 



VttCcY 0.030 e- 2649i 0.048 e a930i 0.032 e~ 3m5i , (31) 
0.848 e 2 ' 339 ' 0.799 e" L353i 0.918 e L838i y 

m WDl) = q Q045 Ge y ^ m {QUDX) = q 10 g Ge y ^ m WDl) = 2 gg Ge y _ (32) 

In Eqs. (|29|). (|31|) the mass matrices are written in a form that better shows deviations 
from the democratic structure. In Eq. (|33j) we give the expression for the CKM matrix, in 
Eq. flMj) the values for the CP parameters p and fj. 
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n 



■ 



8 





q(D 


■ 


q(2) 


■ 


q(3) 




(i) 




u< 2 > 




u< 3 > 




D (l) 




D (2) 




D (3) 



FIG. 1: Profiles of the absolute value of normalized wave functions |^q„|, \ x \j n \ an d |a^£)nl m ^ ne 
theory space (N=10) for the case with Q, U and D localized at n=l. \x^Q n \ with a value of a <C 1 
is less localized. 



V (QUD1) 
V CKM 



0.975 - 0.009 i -0.151- 0160 i -0.001 - 0.003 i 



0.015 + 0.219i -0.669 + 0.709i 0.029 + 0.024i 
\ 0.003 - 0.009 i 0.029 - 0.023 i 0.670 + 0.742 i y 

p«^i) = .12, ^ D1 ) = 0.30, 



with p and r\ defined as 



p = Re{v ud v: b v: d v cb )/\v cd v: b \ 2 , 

fj = im{v ud v: b v; d v cb )/\v cd v; b \ 2 . 



(33) 

(34) 

(35) 
(36) 



• Case 2: Q and U localized in n = 1, D localized in n = N. 





( 


0.918 


e 0.039i 


0.609 e 


-0.590 i 


0.924 e ai35i 


\ 


M^ QU1DN ^ — 66.6 GeV 




0.941 


e 0.038i 


0.637e- 


-0.601 i 


0.946 e°- 132i 






{ 


0.930 


e 0.058i 


0.622 e" 


-0.585 i 


0.935 e ai54i 


/ 



m (QUlDN) = q 002Q ^ m (QUlDN) = q 6g7 ) m (Q^l^) = 16 g 3 Ge y ^ (gg) 
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m a (GeV) m s (GeV) m b (GeV) 



FIG. 2: Solutions for the 6 quark masses in the case with Q, U and D localized at n=l. The 
masses in GeV are evaluated at the Mz scale. The range for each mass is given by the edges of 
the corresponding window. 



Mf U1DN) = 23.2 GeV 



( 0.041 e L959i 0.045 e '°- 025i 0.043 e 2 - 526i \ 
0.037 e L854i 0.042 e°- 003i 0.042 e 2 - 570i 
0.038 e L982i 0.043 e 0M7i 0.043 e 2m5i 



mf U1DN) = 0.0045 GeV, mf ulDN ^ 



0.084 GeV, ml 



(QUIDN) 



2.90 GeV , 



(39) 
(40) 



V 



(QUIDN) 



CRM 



0.975 + 0.029i -0.097 + 0.197 i 0.001 + 0.003i 



; \ 



-0.168- 0.141 i -0.880 + 0A20i 0.039 - 0.011 i 
\ 0.003 + 0.010i 0.039 - 0.007i 0.999 + 0.011 i y 

p(Q^) = .19, r^ 1DJV ) = 0.33. 



(41) 



(42) 



Case 5: Q, U and D localized in n = N. 
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FIG. 3: Solutions for the absolute values of the CKM matrix elements in the case with Q, U and D 
localized at n=l. The range for each element is given by the edges of the corresponding window. 



0.38 - 



0.36 - 



0.34 - 
0.32 - 




0.28 - 



0.26 - 



0.1 0.11 0.12 0.13 0.1 4 0.15 0.16 0.17 0.18 0.19 0.2 
P 



FIG. 4: Solutions for p and r\ in the case with Q, U and D localized at n=l. 
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M (QUDN) = 78AGeV 



( 0.887 e - 4941 0.881 e 0A7Sl 0.913 e a57?l ^ 
0.038 e~ 2 - omi 0.038 e- 2 - 070i 0.041 e~ 2m9i 
0.895 e~ 0AWi 0.877 e-°- 4291 0.929 e~ a316i 
m (QUDN) = 0022 GeV ^ m (QUDN) = 6?4 GeV ? m p^) = i72. 6 Ge \/ , 



/ 



M^ UDN) = 1.37 GeV 



( 0.895 e - 1M9i 0.776 e L668i 0.943 e^ 1 - 622 * \ 
0.058 e 2 - 456i 0.040 e~ L459i 0.063 e 2 - 339i 
0.835 e" 2 - 473i 0.832 e°- 871i 0.893 e~ 2A68i 



m (QUDN) = Q 004g ^ m 



(QUDN) 
s 



0.106 GeV , m[ QraA ° = 2.90 GeV , 



' 0.974 + 0.042 i -0.046 + 0.220 i 0.003 - 0.003 i ^ 
Vckm^ = 0.134- 0.180? -0.676 - 0.701 i 0.020 - 0.033 i 
K -0.010 + 0.006 i 0.022 + 0.030i 0.646 - 0.762? J 

p(QUDN) = 31j -(QUDN) _ g 3Q 
• Case 6: D localized in n — 1, Q and C/ localized in n = N. 



( 0.675 e - L829i 0.824 e" 



Mi D1 ^) = 71.2 GeF 



0.198? 



0.837e 



-0.732 i \ 



V 



0.706 e~ L783i 0.856 e~ ai73i 0.868 e" - 698 ' 
0.671 e- L837i 0.822 e-°- 207i 0.834 e" - 741 ' 
m (DiQUN) = 0Q26 Gey ? m (mQiw) = . 7 25 GeV , m pw) = i 6 9.2 GeV , 



/ 



/ 0.026e°" 6i 0.039e- 3019i 0.044 e - L678 A 
0.027e°- 868i 0.037 e- 3 - 017i 0.042 e - L655i 
0.025 e a957i 0.039 e~ 3 050i 0.044 e" 1 - 710 * 
m (DiQUN) = Q 0Q44 ^ m (Di Q UN) = 088 Gey ? m P W) = 2.91 GeV , 



/ 



tADIQUN) 
V CKM 



-0.972 - 0.075?' -0.069 + 0.213? 0.001 + 0.004?' 



-0.050 - 0.218? 0.974 + 0.013?' -0.039 - 0.016? 
y -0.004- 0.012 i 0.038 -0.014? 0.998 + 0.044?' J 
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We are now ready for comments on the presented numerical solutions. 



(54) 



IV. CONCLUDING COMMENTS 

In this paper we have reconstructed the observed complex mixing of quark flavors, starting 
with the product group Yln=A SU ( 2 ) x U(l)} n at a higher energy scale. The deconstruction 
of this product group into the electroweak gauge group can indeed provide all necessary 
components to generate such mixing. 

We have built a specific models with 20 parameters to fit the quark mass spectrum and the 
CP phase. However, the numerical fit is found only for the "preferred" configurations where 
fermion fields Q and U are localized at the same position in the theory space. Arguably, this 
is because the ratio kjj/k,d of Yukawa couplings can be responsible only for the difference 
in the overall scale of up and down-quark masses, while the more hierarchical internal mass 
spectrum of the up-quark sector (compared to that of the down-quark sector) would still 
require a higher degree of overlapping. 

As far as the structure of mass matrices is concerned, the deviation from democracy is 
moderate. In all the cases, the mass matrices assume a hierarchy with two rows (or two 
columns) having similar absolute value matrix elements, with the third row (or third column) 
having different values, but still similar along that row (or that column). A quite close mass 
matrix structure was found in [18[ , but in a different approach. 

We did not perform a study of the dependence on the number of deconstruction subgroups 
N. We expect anyway that the fitting would be more feasible for larger iV as the wave 
functions and their overlaps then can be tuned more smoothly. In the other direction, the 
constraint from flavor changing neutral current that sets an upper limit on the length of 
extra dimension in the split fermion scenario (see e.g. (2^]) is also expected to set an upper 
limit on the ratio N/V (between N and the VEV of link field) in the deconstruction theory. 
We however leave a more careful analysis of these and other relevant phenomenological issues 
for future publications. 
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APPENDIX A: FERMION ZERO MODE IN DIMENSIONAL DECONSTRUC- 
TION 

In this appendix we will work out the general expression of zero eigenstate of the matrix 
of the type (JBJ). This mode plays a special role because it will be identified with the SM 
chiral fermions. To simplify the writing, here we denote this zero eigenstate generally as 
{xi,x 2 , • • • , %n-i} while in Section III we will restore all omitted scripts Q, U, D, 

1. Zero-mode localization at the end-point n=l 

The equation set determining the zero eigenstate (0) is 



V\ 2 xi - MVx 2 + V 2 x 3 = <£> V*x x - Mx 2 + Vx 3 = 
-MV*x 1 + (M 2 + \V\ 2 )x 2 - 2MVx 3 + V 2 x 4 = 



(Al) 



(A2) 



V* 2 x l - 2MV*x 2 + (M 2 + 2|Vf )x 3 - 2MVx A + V 2 x 5 = 



(A3) 



V* 2 x N ^ - 2M\"r x , + (M 2 + 2|Vf )xiv-3 - 2MVx N „ 2 + V 2 x N . 1 = 
V* 2 x N . 4 - 2MV*x N ^ 3 + (M 2 + 2|Vf )xiv-2 + (V 2 - 2MV)x N ^ = 
V* 2 x N - 3 (V* 2 - 2MV*)x N _ 2 + [M 2 - M(V + V*) + 2\V\ 2 ]x N . 1 = . 



(A4) 



(A5) 



(A6) 



After a bit of algebra, we can equivalently transform this equation set into 




(A7) 



(A8) 



Xn-3 — X N _ 2 — \p\ X N -i 



(A9) 



18 



X 



N-2 — 



pX 



N-l i 



(AlO) 



where we have introduced new parameter and variables 

v \v\e 



■w 



P = \P\e = 

X n ee (p*)^-"- 1 



(All) 
(A12) 



M M 

(n = l,..,N-l). 

We note that V (and p) is a complex parameter in general (see appendix B). The new simple 
recursion relation allows us to analytically determine the set {Xx, ..,X^-i} (and then the 
zero eigenstateixi, ..,xn-i} ) for any p (i.e. for any real M and complex V). After some 
combinatorics l25|] we obtain for 1 < n < N — 3 



X n — 

^(N-3-n-k\ . , —(N-2-n-k. 

E (-H 2 ) fc+1 ^-i + E l(-H 2 )^- 2 

k=o \ k k=o \ k 



(A13) 



(N-3-n-k)\ 2k+1 
' £^ Q k\{N -3-n-2k)V lPl ' 

and for n = N - 2 (see (jXTHjl ) 



(N — 2 — n — k)\ 
,^k\{N-2-n-2k)\ 



Xn~i + E Tim — n — i — oT\t( — \p\ 2 ) k XN-j ■ 



Xn-2 — (1 — p)ATjv_i . 

Using the equality 

m \ I m \ I m + 1 
p I \p+l J \ p+1 
we can rewrite (|A13|) as (with 1 < n < N — 3) 



(A14) 



(A15) 



X„ 



E 

fc=0 



N — 1 — n — k 
k 



IpIY-pE 

fc=0 



N -2-n-k 



2\k 



HpH 



X N -l (A16) 



Again, using another equality 



24] 



sinhpx = sinhx E(~V 



fc=0 



p — 1 — k 



(2 cosh x 



p-l-2fe 



(A17) 



V 



we obtain for 1 < n < N — 3 and for \p\ < | 
2Xtv_! 



A. 



1 -4|p| 



p\ N ~ n sinh (A" - n)a - p|pr sinh (X - 1 - n)a 



(A18) 
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with 



cosh a = —r-r a = cosh 1 



(a > 0) 



(A19) 



2|p| 2\p\ 

For |p| > ~, the expression of Xn is similar to (jA18|) but with hyperbolic functions (sinh 
and cosh) replaced respectively by trigonometric ones (sin and cos). 
Finally, from (jA12|) we have altogether 



-ind 



sinh (N - n)a - e ie sinh (N - 1 - n)a (1 < n < N - 1) , (A20) 



where C is the normalization constant determined by the normalization equation 



N-l 

E 

71=1 



1 . 



(A21) 



We note that this normalization is nothing other than the unitarity condition of the rotation 
matrix U (see (IH|) and (l23|) . (12^ ). 



2. Zero-mode localization at the end-point n=N 

The chiral boundary conditions (CBC) and the value of parameter \p\ = \V\/\M\ are two 
crucial factors that determine the localization pattern of the chiral zero-mode of fermion. 
For e.g. in the previous subsection we have seen that, when |p| < 1/2, along with CBCs 
Qir — Qnr = , 07V-1 nQn,l = VqQn-i,l @ we can localize the left-handed zero mode 
of Q field around site n = 1 (lA"20l) . 

On the intuitive ground, we expect that the "mirror image" of (J2J) (apart from the re- 
quirement \p\ < 1/2) 

Qi,r = Qn,r = , <t> Q { 2 Qi, L = VqQ2,l , (A22) 

would produce a left-handed zero mode of Q field localized at n = N. A similar calculation 
indeed confirms this localization pattern. Specifically, if we denote t/i (i = 2, .., N) the zero- 
mode subject to CBCs (jA22|) . and Xj (j = 1, ..,N — 1) subject to CBCs (J2J) as before, we 
find 

y i = x* N+1 _ i (i = 2,..,N), (A23) 
or even more explicitly (see ()A20|) ) 

Vn = Ce i{N+1 - n)e [sinh (n - l)a - e~ ie sinh {n - 2) a] (2 < n < N) . (A24) 
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APPENDIX B: THE COMPLEX- VALUED LINK FIELD VEV FROM BROKEN 
WILSON LINE 



Since the link fields n>n +i transform non-trivially under two different groups, we may 
expect its VEV to be complex in general. It is because in this case the VEV's phase could 
not be rotated away in general. The standard and rigorous method to determine the VEV 



2J that there 



is to write down and then minimize the corresponding potential. It turns out 
always exist ranges of potential parameters which generate complex VEV. In this appendix, 
however, we just recapitulate the complexity nature of link field VEV from the latticized 
extra dimension perception which is derived in in details. Though the approach taken 
in this work does not strictly stem from latticizing the fifth dimension, this perception could 
still serve as the principle illustration. 

To make the connection between DD theory and its latticized ED counterpart, we inter- 
pret the link field as a Wilson line connecting two neighboring branes 

r(n+l)a 

<t>n,n+i ~ exp( / igXydy) ~ exp(^axn) , (Bl) 

J na 

where Xn essentially is the ED component of gauge field, g and a are gauge coupling and 
lattice spacing respectively. 

Following the DD symmetry breaking [SU(2) x U(1)} N — > [SU(2) x U(l)}, only one linear 
combination xo of link fields remains massless at the classical level 

1 N 

Xo = ^=Y,Xn. (B2) 
V n=1 

In the leading order with radiative correction, by minimizing the 1-loop effective potential 
of xo, one obtains 

(Xo) = (keN). (B3) 

gay N 

From (jBljl , ()B2j) , one see that in the leading order the link fields assume a uniform complex 
VEV 

, . j 2 2 JvTT , , . 

(<Pn,n+l) ~ eX P(^r) • ( B4 ) 

Actually, this phase can be considered arbitrary. 
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APPENDIX C: WAVE FUNCTION OVERLAP IN THEORY SPACE 



In this appendix we present the analytical expressions of zero-mode wave function overlaps 
in the theory space, from which follow the mass matrix elements M^f (J27)l . (|2Sjl . These 
expressions in turn were compiled using the exact solutions (|A20j) . (jA24|) for the wave 
functions. In what follows we use XX to denote the overlap of two wave functions localized 
at the same site n=l, and XY the overlap of the first wave function localized at n=l and 
the second at n=N. All other overlap configurations can be easily found by virtue of relation 
(1X23J) . 

It follows from Eqs. JZ@, (EHD that 



xx = E« 



\ n=l 



CCr. ( P (N-l)(ie 1 -ie 2 -a 1 -a 2 ) _ 1 

4 I l _ e -(i9 1 -i0 2 -ui-a 2 ) ^ ^ ' 

e (AT-l)(i6»i-i9 2 -Q!i+a2) _ I 



1 


_ g— (i6\— i0 2 — 01+02) 




e (N- 


-l)(idi—i9 2 +ai-a 2 ) 


1 


1 


_ g— (261—46*2+01-02) 


— 




-l){iei-i$2+ai+a 2 ) _ 


1 



N(a 1 -a 2 ) /-^ _ g-iei-ai^-^ _ ^id 2 +a 2 ^ 



e N(-ai+a 2 ) /-^ _ g~ i6, i+ Q i)(;L _ g i6, 2-a!2' 



I p -N(ai+a 2 )(-i _ -i8 1 +a 1 \(-t _ i0 2 +a 2 \ 1 

" t ~ 1 _ g-(^i-ie 2 +ai+a 2 ) e ^ n x e n 

+ C x C 2 e< N - l ^- e ^ sinh «! sinh a 2 , (CI) 

JV-1 \ 

(1)* „,(2) 



/JV-1 \ 

YY = ftf+feW 

\n=2 ) 



x n-iVn 



= C^e 101 - 292 [sinh (Y - l)ai - e"^ 1 sinh (Y - 2)ai] sinh a 2 

(7,f7, / g(Af-2)(iei-i6» 2 -ai-a 2 ) _ 1 

4 I 1 — e -(*%-*02-«i-a2) ^ ^ ' 

(Ar-2)(i6»i-i02-«i+a 2 ) _ 1 

I _ g-(jei-i9 2 -ai+02) ^ ' 

(Af-2)(jei-i6»2+ai-a 2 ) _ 1 
1 e i9 1 -i6 2 +a 1 -a 2 e -Nax+a 2 _ e ~i9 1 +a 1 \ _ g-r02+o 2 N 

1 g— («0i— 16*2+01-02) \ /\ / 

(Ar-2)(i6>i-«0 2 +ai+a2) _ 1 \ 
]_ _ g-(iei-i6 l 2+ai+02) ^ '] 

+ C^e^- 1 ^ 1 -™ 2 sinh ai [sinh (Y - l)a 2 - e^ 2 sinh (Y - 2)a 2 ] • (C2) 

where Ci,C 2 are the normalization factors, which are determined also from the overlap of 
the respective wave function with itself. 
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APPENDIX D: NUMERICAL TABLES 

TABLE I: Central values and uncertainties for the masses of the 6 quarks evaluated at Mz, for the 
two ratios m u /m d and m s /m d , for the absolute values of the CKM matrix elements and the CP 
parameters p, fj 



x . <Xi> \x™ ax - x™ in \/2 



m u 2.33 x 1(T 3 0.45 x 1(T 3 

m c 0.685 0.061 

m t 181 13 

m d 4.69 x 10~ 3 0.66 x 10~ 3 

m s 0.0934 0.0130 

m b 3.00 0.11 

m u /m d 0.497 0.119 

m s /m d 19.9 3.9 

\V ud \ 0.97485 0.00075 

\V US \ 0.2225 0.0035 

\V ub \ 0.00365 0.0115 

\V cd \ 0.2225 0.0035 

\V CS \ 0.9740 0.0008 

\V cb \ 0.041 0.003 

\Vtd\ 0.009 0.005 

\V U \ 0.0405 0.0035 

|Vtb| 0.99915 0.00015 

p 0.22 0.10 

fj 0.35 0.05 
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TABLE II: 20-parameter space solutions found in 4 different cases of the model presented in Sec. 
IIIA (N=10 for all cases) 





(QUD1) 


(QU1DN) 


(QUDN) 


(D1QUN) 




2.290 


0.208 


2.311 


0.215 


«Q2 


0.007 


0.236 


0.011 


0.251 


«Q3 


1.497 


0.220 


1.656 


0.213 


am 


0.771 


0.439 


0.623 


0.572 




0.759 


1.900 


0.603 


0.823 


"C73 


0.927 


0.433 


0.722 


0.745 


am 


0.829 


0.027 


0.794 


0.541 


OID2 


0.535 


0.022 


0.471 


0.057 


&D3 


1.123 


0.029 


1.105 


0.022 


Oqi 


1.002 


0.974 


1.234 


0.959 


Q 2 


3.862 


0.983 


4.702 


0.958 




1.678 


0.989 


1.999 


0.965 


Ojji 


1.163 


1.407 


1.276 


0.197 


dm 


1.100 


14.78 


1.246 


13.67 


&U3 


1.247 


-5.290 


1.420 


-5.52 


Odi 


-0.204 


-0.535 


-0.135 


-0.182 


&D2 


3.098 


4.715 


2.809 


3.471 


&D3 


0.086 


9.645 


-0.174 


9.499 


kuv/V~2 


78.37 


66.63 


78.36 


71.21 


KDV/V2 


1.35 


23.24 


1.37 


26.48 



[1] N. Arkani-Hamed, A. G. Cohen and H. Georgi, Phys. Rev. Lett. 86, 4757 (2001). 
[2] C. T. Hill, S. Pokorski and J. Wang, Phys. Rev. D 64, 105005 (2001). 



24 



[3] N. Jackiw and S. Rebbi, Phys. Rev. D 13, 3398 (1976). 

[4] N. Arkani-Hamed and M. Schmaltz, Phys. Rev. D 61, 033005 (2000) 

[5] P. Q. Hung and N-K. Tran, Phys. Rev. D 69, 064003 (2004). 

[6] N. Cosme, J. M. Frere and L. L. Honorez, Phys. Rev. D 68, 096001 (2003). 

[7] N. Cosme and J. M. Frere, Phys. Rev. D 69, 036003 (2004). 



[8] B. Gdzadkowski and J. Wudka, |hep-ph/0401232 



[9] Y. Hosotani, Phys. Lett. B 129, 193 (1983). 
[10] E. A. Mirabelli and M. Schmaltz, Phys. Rev. D 61, 113011 (2000). 
[11] G. C. Branco, A. De Gouvea and M. N. Rebelo, Phys. Lett. B 506, 115 (2001). 
[12] A. Soddu ans N-K. Tran, Phys. Rev. D 69,015010 (2004). 
[13] W. Skiba and D. Smith, Phys. Rev. D 65, 095002 (2002). 
[14] C. T. Hill and A. Leibovich, Phys. Rev. D 66, 016006 (2002). 
[15] H. Abe, T. Kobayashi, N. Maru and K. Yoshioka Phys. Rev. D 67, 045019 (2003). 
[16] P. Q. Hung and N-K. Tran, Phys. Rev. D 70, 126014 (2004). 
[17] P. Q. Hung and M. Seco, Nucl. Phys. B 653, 123 (2003). 
[18] P. Q. Hung, M. Seco and A. Soddu, Nucl. Phys. B 692, 83 (2004). 
[19] S. Nojiri, S. D. Odintsov and A. Sugamoto, Phys. Lett. B 590, 239 (2004). 
[20] S. Kirkpatrick, C. D. Gelatt and M. P. Vecchi, Science 220, 671 (1983); S. Kirkpatrick, J. 

Stat. Phys. 34, 975 (1984); N. Metropolis, A. Rosenbluth. M. Rosenbluth, A. Teller and E. 

Teller, J. Chem. Phys. 21, 1087 (1953). 
[21] W. H. Press, S. A. Teukolosky, W. T. Vetterling and B. P. Flannery, Numerical Recipes in C, 

2nd ed. (Cambridge University Press, New York, 1992). 
[22] F. Bauer, M. Lindner and G. Seidl, J. High Energy Phys. 0405, 026 (2004). 
[23] A. Delgado, A. Pomarol and M. Quiros, J. High Energy Phys. 0001, 030 (2000). 
[24] I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series and Products, 5th ed. (Academic 

Press, Boston, 1994). 

[25] Another somewhat simpler solution of the above equation set is presented in ^(| 

[26] The finiteness of 1-loop effective potential requires the mass of fermionic tower be trigonometric 



function of the mode number 



14]. In the continuum limit such as in a "S l /Z2 11 compactifica- 



tion, the orbifold boundary conditions (J2J) can fulfill this requirement. 



25 



